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Abstract. We perform a calculation of the linewidth of a micromaser, using the master 
equation and the quantum regression approach. A 'dephasing' contribution is identified 
from pumping processes that conserve the photon number and do not appear in the photon 
statistics. We work out examples for a single-atom maser with a precisely controlled cou- 
pling and for a laser where the interaction time is broadly distributed. In the latter case, we 
also assess the convergence of a recently developed uniform Lindblad approximation to the 
master equation; it is relatively slow. 

1 Introduction 

The linewidth of a laser and the spectrum of a resonantly driven two-level system are paradigmatic ex- 
amples that probably spawned the field of quantum optics |[Tl|2l[3l. In both cases, one has to deal with 
an open quantum system that is driven towards a steady state by the balance between pumping and dis- 
sipation. The characteristic features are intimately linked to the quantized nature of the electromagnetic 
field, with spontaneous processes determining the (laser) linewidth and the triplet spectrum (resonance 
fluorescence). 

A fertile playground to investigate these issues in great detail is the micromaser or single-atom 
laser EUs). Its basic ingredients are well known: a high-Q cavity sustaining an electromagnetic field 
mode, and a dilute atomic beam that delivers energy by entering the cavity in an electronically excited 
state. In this device, a delicate control of the interaction time is possible and permits the observation of 
non-classical behaviour of the cavity field like quantum jumps ||6l or the coUapse into photon number 
states i7"8l. In addition, the atoms leaving the cavity provide information that can be used, for example, 
to measure in a non-destructive way the cavity photon number |9, .10J . 

In this paper, we discuss the linewidth of the micromaser or in other words, the coherence time of 
the cavity field. This has been studied a few years after the introduction of the basic principles of the mi- 
cromaser 1 12J by Scully and co-workers [13J . The basic ideas are simple: the maser (or laser) spectrum 
is defined by the Fourier transform of the two-time correlation (or coherence) function of the cavity 
field operator 114). To calculate these correlations, one invokes the so-called 'quantum regression for- 
mula' II15II16I that provides a link to the time evolution of the system density operator. This evolution 
is conveniently given in terms of a master equation whose eigenstates and eigenvalues determine the 
maser spectrum. Their computation typically requires numerical work, and it is difficult to find a physi- 
cal interpretation of the results. In particular, more than one eigenvalue contributes to the result around 
the trapping states where phase transitions happen in the micromaser |17). We suggest here a formula- 
tion of the maser linewidth (i) that provides a simple physical interpretation; (ii) that can be computed 
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fairly easily once the steady state solution to the master equation is known; and (iii) that interpolates in 
a natural way between different eigenvalues of the master equation. Our methods are similar to those 
used by Scully 1. 1 3] and McGowan and Schieve [18J where a iinewidth operator' that depends on the 
photon number is identified, whose expectation value with respect to the steady-state density operator 
yields the Iinewidth. The main difference of our approach is its straightforward formulation and sim- 
plicity. The results are generally in good agreement with previous work. We predict a weak oscillation 
of the Iinewidth as a function of the pumping strength that is related to a dephasing process related to 
the non-destructive monitoring of the cavity photon number by atoms that do not deposit their excitation 
energy in the cavity. 

The outline is as follows. We fix the notation in Sec 12. II introduce our formula for the Iinewidth 
(Eq.Q of Sec l2.2l ) and work out the details for two cases: a 'perfect maser' where the interaction time r 
is fixed (Sec l3.1l l and a 'conventional laser' where r is exponentially distributed (Sec 13. 2b . We conclude 
by a comparison to a recently developed approximation [191 that turns out to describe the laser above 
threshold fairly well, but becomes increasingly complex in the strong pumping regime (Sec l3.3b . The 
Appendix contains some technical details of the calculation of the micromaser Iinewidth. 

2 The micromaser and its Iinewidth 
2.1 Master equation 

The master equation for the micromaser we use here is of the form ll20l l3l 



The first two lines describe cavity damping and its thermal excitation by the surrounding in the usual 
way. The thermal occupation number nth depends on temperature and gives the average photon number 
at equilibrium, without pumping. The maser is pumped by a stream of excited two-level atoms that enter 
the cavity one by one with Poissonian statistics (rate r). The atoms interact with the maser mode for 
a time r that is distributed according to the measure dp(T), and g is the one-photon Rabi frequency 
for this interaction. The operator ip = (aat)^/^ and sinc(x) = sin(2;)/a;. The conventional pumping 
parameter is given by 6* = {r/ ny/'^gr. 

In the strong coupling regime, non-classical states of the radiation field like number states are gener- 
ated whenever the photon number n satisfies grV n + 1 = tt, 27r, . . . lil2J . Then, the last term in Eq.([T]i 
is zero and the photon number in the cavity remains unchanged. Physically speaking, each incoming 
atom performs an integer number of full Rabi cycles and leaves the cavity in the excited state again. 
This perfect scenario is perturbed by the fluctuations of the interaction time r and a not perfectly con- 
trolled coupling strength g that depends, for example, on the atom crossing the laser mode at a node or 
an anti-node. A typical choice for the distribution dp{T) in conventional lasers is an exponential one, its 
mean value f giving the lifetime of the excited atomic states. 

2.2 Maser spectrum 

The maser spectrum is defined by the Fourier transform of the two-time correlation function il3l [T4l 



dp 
dt 




(1) 
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and is usually computed in the stationary regime where the correlation function only depends on the 
time difference t with respect to which the Fourier transform is taken. We therefore drop the to argument 
from Eq.(|2]l. Alternative, but equivalent definitions are based on the decay of the average electric field 
operator {a{t)) (Ref.| 13|) or the correlation function of e"^'*-' with (/) being the phase operator (Re f. 1211 ). 

Using arguments in the spirit of the quantum regression formula, one shows that the coiTelation 
function can be computed as the expectation value B14ll21lfr7ll 

g{t) = tr [a^P{t)] (3) 

of a photon creation operator with respect to a 'skew density operator' P. The latter obeys an equation 
of motion identical to the master equation ([U, but with the 'initial condition' P(0) = apss where pss is 
the stationary solution that turns out to be diagonal in the photon number basis. A trivial consequence 
is the normalization to the average photon number g(0) = tr (npss) = {fi)ss- 

The expansion of P{t) in eigenvectors of the master equation yields a sum of exponentials in time 
and hence of Lorentzian spectra. Usually, the spectra are computed in the frequency domain and their 
characteristic width computed numerically. We adopt here a different approach that gives straightfor- 
ward analytical results, but restricted to the regime where at most two eigenvalues are dominating. We 
emphasize this happens over a relatively large part of the parameter space (see Ref. lJTJ ). Our definition 
of the maser linewidth D is 

D 1 dg 



2 (n).s dt 



(4) 





where the time derivative is computed directly from the master equation, using the quantum regression 
formula (O. To justify this choice, note that the eigenvector expansion for the correlation function is of 
the form 



3 

where the prefactors satisfy the sum rule J^j 9j — 5(0) = (^)ss- This leads to 

D = (6) 

hence an arithmetic mean of the exact eigenvalues. Now, it appears from Ref. llTTl that the relative 
weights gj/ {n)ss are often close to either zero or one. Only when the two lowest eigenvalues are cross- 
ing, both are contributing significantly. If a single eigenvalue is relevant, our linewidth D coincides with 
it; if two contribute, Eq.® is just another convenient way to combine the two into a single parameter. 
We note that the conventional full width at half maximum (FWHM) for a sum of two Lorentzians of 
equal weights, but of different widths, leads to the geometric mean of the widths. This significantly dif- 
fers from the arithmetic mean only when the widths are very different. But, as mentioned above, equal 
weights occur typically when the eigenvalues cross. 



2.3 Explicit formula and comparison 

The time derivative in Eq.(|4]l leads to, using the regression formula Q 



D 



-tr 



dP 
'dt 



(7) 



This expression involves the master equation for the skew operator P which can be calculated easily by 
adopting a Lindblad form. One gets the familiar result 



D(n)^^ = - 



A I aPss + 



4, at 



Lxapss^ 



(8) 
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where the L\ axe. the so-called jump or Lindblad operators. For the micromaser, they can be written as 
(the index is running over A = 0, 1, (c, k), (s, k)) 



Lq = \J k(1 + nth) a (9) 

il = ^JK.nt\,a} (10) 

ic,fe = \/rAp(Tk) cos(gTfc^) + 1/^ (11) 

LsM = \/r Ap{Tk) gTk ahmc{gTk(p) (12) 

where we have broken the integral over r into a Riemann sum with weights Ap{Tk). The term pro- 
portional to the unit operator in Lc.k actually does not contribute to the Lindblad master equation; an 
expression for the real number fk can be found in Ref.| 19|, it is not needed here. (Note that the term 
—rp in Eq.O is generated by the anticommutators in the Lindblad form, after performing the Riemann 
integral.) 

We immediately see that one gets the hnewidth as the average value of some operators with respect 
to the steady-state photon statistics, similar to what has been found previously M13II18I . In Appendix]?] 
these operators are computed and simplified using the detailed balance condition for the steady state. 
The result is 

D{h)ss = unth + ^ Ap{Tk)tT |r( cos{gTk(p) - cos{gTkh^^^))^hpss 

k 

+ r{ipsm{gTk(p) - n^^^ sm{gTkn^^^)y p^sj (13) 



where h = a'' a is the photon number operator and we have kept for the ease of comparison the discrete 
summation. 

This result should be compared to previous calculations. In their seminal paper llT3l . Scully and 
co-workers make the approximation of a sufficiently narrow distribution of photon numbers and get a 
linewidth 

^ _ l_f2nth 2 I — — ] (14) 

McGowan and Schieve improve this calculation by allowing for multiple eigenvalues of the Liouville 
operator |18|. They find an expression that gives a good approximation (to a few percent) to the ex- 
act eigenvalue calculation. The method is using a numerical evaluation of the FWHM of the (non- 
Lorentzian) spectrum which we want to avoid here. But their approximation can also be formulated in 
terms of the correlation function g{t), and applying our definition Q, one has 

D{h)ss ~ 2 ^ Ap{Tk)ti (r -r sm(gTk(p) sm{gTkn^/'^) - r cos{gTk(p) cos{gTkn^/'^) 



+ k(1 + nth)(n - 5 - [n{n - 1)]^''^) + Knth{n + 5 - (ph^^^))hps, 
where we have added the average over r as before. 



(15) 



3 Discussion 
3.1 Perfect maser 

To get a better understanding of the physics behind the maser linewidth, we first fix the interaction time 
(no fluctuations) and plot in Fig[T|the contribution of the Fock state |n) to the linewidth, that can be read 
off from the operators under the trace in Eqs. (|13|15l ). In the same plot, we also give the photon statistics 
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Fig. 1. Contribution to the maser linewidth for a given photon number. We plot the n-dependent term under the trace 
of Eg.fTst (undulating curve) and Eg.JlSt ('McGowan') . The normalization is arbitrary, and the interaction time is 
fixed (no average). The peaked curves give the steady-state photon statistics. Solid and dashed curves correspond 
to two pumping parameters: 9 = 2.1 n (gr ~ 0.47, left peak, solid line) and 6 — 2.2 n {gr ~ 0.49, right peak, 
dashed line). We fix the pumping rate to r = 200 k and nth = 0.1. 



(peaked curves). The chosen parameters correspond to the first trapping transition just above 6* w 27r 
where the maximum of the photon statistics jumps from lower to higher values. We see that Eqs.lfTSll 
predicts a weak oscillation of the linewidth as the mean photon number increases, while Eq.dTSb is 
essentially flat in this regime. These oscillations arise from the different prefactors (p and fi^/"^ of the 
sine terms in Eq.lfTsl). Indeed, performing and expansion in 1/n, the trigonometric functions can be 
simplified and give an expression identical in leading order to Eq.lfTSll. The remainder is of relative over 
\/tt? and contains the sinusoidal oscillation seen in Fig[T] 

In Fig|2l we plot the linewidth as a function of the pump parameter 6 and see that all three expres- 
sions give very similar values for D. We also note (see inset) that the product D{fi)ss/ n that measures 
the deviation from the Schawlow-Townes limit of a conventional laser is essentially a quadratic function 
of 9, similar to the formula found in Ref.|,13 |: 

^ k9^ + k{1 + 2nth) ,,,, 
D{n}ss ~ (16) 

It is interesting to note that the sharp features of D near the trapping states essentially arise from the 
corresponding oscillations in the average photon number. 



3.2 Exponentially distributed interaction time 

Our formula (fT3] ) can be averaged explicitly over an exponentially distributed interaction time r. The 
result takes the form (for simplicity, we make the replacement gf i-^ g) 

_ , [l+g^(3n + 2)+g^(n + l)(4n+l)] p^s 
^'"'^^ ~ [1 + mh + 1)][1 + 2g^2h + 1) + g4] '> 

The corresponding photon statistics is given by the recurrence relation 

2P 

1 + 4g^(n + 1) 

where 9 = 5f(r/K)^/^ is the pumping parameter corresponding to f. Here, we have neglected the 
thermal effects and set nth = 0. 

The linewidth ( fTTI l is plotted in Figl3] as a function of the (average) pump parameter 9. The sharp 
features at the trapping transitions are smoothed due to the averaging over r. Above the threshold. 
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Fig. 2. Comparing expressions for the micromaser linewidth vs. the pumping parameter 9 = (r j k)^^"^ gr . No 
average over the interaction time r here. Main plot: the expressions by Scully rEq.l ll4t . dashed line] and by Mc- 
Gowan and Schieve [Eg.lllSt, red triangles] practically coincide, except far below threshold (S/tt <^ I/tt). Eq.J13l> 
suggested here predicts in addition weak oscillations. Inset: linewidth normalized to the Schawlow-Townes value, 
D{n)ss/K vs. pumping parameter Q around the first trapping state. A similar smooth behaviour is found also at 
larger values of Q. 

Parameters: r = 50 nth = 0.01. We scan Q via the product gr. 




9 >1/ -\/2, the linewidth is quadratic in 9. This can be derived from the Scully formula ( fT4b by assuming 
that the average photon number is large and expanding the sine to lowest order This gives a term 
''^{qt)'^ I ^{f>)sB = K9'^/4{n)ss, and taking the average over an exponential distribution for r, we get 

^-ifr^ (.9) 

This is plotted in dashed in FigUand gives good agreement way above threshold. 



3.3 Uniform approximation to the master equation 

We also give in Figl3]the results from a recently developed approximation to the laser master equation 
that is built to be valid beyond the weak coupling regime gf ^ 1, see Ref.|19|. The basic idea is to 
replace the large set of Lindblad operators in Eqs.(l9]l-(fT2]i by a smaller one, using an expansion adapted 
to the probability measure dp(r). More specifically, we represent the operators cos{gTip) and sin(gri^) 
as sums of normalized orthogonal polynomials in r/f , with coefficients that are rational functions of 
gf(p. The expansion is truncated at polynomials of a given order and gives good agreement for the 
average photon number and its variance. The agreement persists above the laser threshold where a 
simple expansion in gf fails, as illustrated in Ref.[191. Fig|3]shows, however, that for the laser linewidth, 
polynomials of fairly high order (up to the seventh) have to be taken into account to describe the strong 
pumping regime 9 ^ 1. In fact, if we compare the n-dependent linewidth (inside the trace of Eq.(l3]i) 
with its uniform expansion, we see a good agreement only up to photon numbers where gfn^^^ is 
comparable to the maximum retained order. Whenever the photon statistics covers photon numbers 
beyond that value, the (truncated) uniform approximation breaks down. This can be translated into a 
upper limit for the pumping parameter that scales linearly with the order — a behaviour that can be 
qualitatively seen in Figl3] 



4 Conclusions 



We have discussed in this paper an alternative calculation of the linewidth of a laser, with emphasis on 
the micromaser. The expression we suggest can be immediately computed from the master equation. 
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Fig. 3. Linewidth for a laser with exponentially distributed coupling parameter gr. Solid line from Ea.J17t. Dashed 
line: Ea.U9t from Scully and co-workers. Red triangles: uniform approximation (see Ref.[ 19|) up to order 7, 3, 1 
(top to bottom) in gfif> for the expansion of sm{gTi^) (see Sec l3.3t . For the cosine, the maximum order is 6, 2, 0, 
respectively. 

Parameters: gf = 0.3 and 6 — gf{r/n) is scanned via the pumping rate r. The arrow marks the laser threshold 
6 — 1/ \p2 where the small-signal gain 2r{gf)'^ equals the loss rate k (see Ref.t3il). 

using the quantum regression theorem, see Eq.(|3]i. The result emerges as a average over the photon 
number statistics, similar to previous work that we find good agreement with. 

In our expression to the Unewidth, a genuine contribution from a 'dephasing' mechanism is evident. 
It is due to the cosine terms in Eq.(fT3]l that can be traced back to those pumping events where the 
pumping atom is detected again in the excited state after crossing the cavity. These processes cannot 
contribute to the laser gain, this is why only the sine terms appear in the photon statistics [Eq.(l24l)l. 
The combined atom-ncavity system does get affected, however, by a signed probability amplitude. This 
is very similar to non-destructive measurements of the photon number using atoms in superposition 
states [TOl. It is the average over these signed probability amplitudes, since they depend on the photon 
number, that contributes to the maser linewidth. Our formula for the linewidth nicely embodies this 
mechanism via the differences cos{gT\/n + 1) — cos(grY^) that occur in Eq.(fT3]l. Even when the 
interaction time is broadly distributed, dephasing contributes a slight increase of the linewidth near the 
lasing threshold (FigO. 

Acknowledgements. I thank Giovanna Morigi and the colleagues from the quantum optics work- 
shop in Palermo (June 2007) for instructive discussions. 

Appendix 

The contribution of the thermal absorption and emission terms to the linewidth is straightforward using 
the commutators \a) , Lq\ = — + rttti) and [a, Li] = ^k, nth- This gives 



or Z3|j.j^ = K. (Perfectly reasonable: only spontaneous, no stimulated processes contribute to the linewidth.) 
The pumping terms can be simpUfied with the help of the identity 



valid for any smooth function F{(p). For example, commutators are transformed into differences 




(20) 



a 



)F{ip)^F{h^'^)a), 



(21) 




(22) 



8 



Will be inserted by the editor 



We thus get from the cosine terms 

^Icos (")ss = X! ^P('^fe)tr I (cos{gTk(f) - cos(.grfcri^/^)^ hpss \ (23) 

The result from the sine terms is sHghtly more comphcated. We now use the detailed balance for the 
photon statistics to bring it into a manifestly positive form. Recall that the steady state density matrix 
has diagonal elements p„ = {n\pss\n) that satisfy the recurrence relation lfT2l 

Pn+i = — {nth + O'^smc'^igrVnTl)) Pn (24) 

Including the average over the interaction time r, we find 

= + 71th)Pn+l + KJ^thPn + K Z\p(Tfc )6>fcSinC^ (ffTfc Vn + l)pn (25) 

fe 

with the pumping parameter 9k = {r/Ky^^{gTk ). We multiply this expression with n + 1, sum over n 
and subtract it from the linewidth. There is partial cancellation of the terms proportional to k with the 
thermal contribution ( l20l i. and finally only the thermal occupation remains, as stated in Eq.(fT3]l. 
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